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In the present letter we suggest a new theoretical model for a quantitative de- 
scription of the magnetoreception mechanism in birds. The considered mechanism 
involves two types of iron minerals (magnetite and maghemite) which were found 
in subcellular compartments within sensory dendrites of the upper beak of several 
bird species. The analysis of forces acting between the iron particles shows that 
the orientation of the external geomagnetic field can significantly change the prob- 
ability of the mechanosensitive ion channels opening and closing. The performed 
theoretical analysis shows that the suggested magnetoreceptor system might be a 
sensitive biological magnetometer providing an essential part of the magnetic map 
for navigation. 

A large variety of animals possess a magnetic sense. The best-studied example is the use 
of the geomagnetic field by migratory birds for orientation and navigation during migration. 



Reviews of these studies are given in Refs. 



2| . In the present letter we address this problem 



from the theoretical point of view. Based on the known experimental observations we 
develop a physical model which we use for the description of magnetoreception phenomena 
in birds. The suggested model is based on the interaction of magnetic particles consisting 



of ferrimagnetic iron-minerals magnetite (Fe^O^ 
observed in the beak of different bird species 



n n n 



and maghemite (7— Fe203), which were 



Based on the analysis of forces which 



act between these particles we show that the considered iron-mineral system can serve as a 
magnetoreceptor with distinct orientational properties. We demonstrate that- depending on 
the orientation of the external magnetic field- the probability of opening of mechanosensitive 
ion channels significantly changes, thus leading to different nerve signals. The nerve signals 
are delivered to the brain causing a certain orientational behavior of the bird. 

The histology studies of the upper beak of homing pigeons 0, 0] showed that iron min- 
erals are concentrated in six symmetrical spots near the lateral margin of the skin of the 



maghemite 
platelets 




FIG. 1: Magnetoreceptor unit consisting of ten maghemite platelets (boxes) and a magnetite 
cluster (sphere). The coordinate frame shown here is used in the computations. The direction 
of the external magnetic induction vector B is characterized by two polar angles <3? and 0. The 
magnetic moments of the maghemite platelet i, rhi, and of the magnetite cluster (M) are indicated 

upper beak inside the dendrites of nerve cells. For the study of the magnetoreception func- 
tion of the dendrite a primary magnetoreceptor unit has been defined, being the smallest 
structure possessing the magnetoreception properties of the whole dendrite. The magne- 
toreceptor unit consists of ten maghemite platelets and one magnetite cluster as shown in 
Fig. [H Experimental observations 0, Q suggest that the dendrite contains about 10-15 
magnetoreceptor units, which should have similar behavior in the external magnetic field. 
Therefore, if the entire dendrite is subject to the external magnetic field the repetition of the 
magnetoreceptory units increases the functional safety of the dendrite magnetoreception. 



The geometry of the magnetoreceptor unit is determined from the experimental obser- 
vations |5j. Thus, the maghemite platelets have the dimensions 1 x 0.1 x 1 /im and the 
magnetite cluster has the diameter of 1 /im. The maghemite platelets are located in the 
(xz)-plane being aligned along the x-axis (see Fig. [T]). The distance between two neighboring 
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platelets is equal to 0.1 jim.. 

The size of a single maghemite platelet (MP) is sufficient for the formation of magnetic 
domains in the (xz)-plane of the platelet |6( (see Fig. [1]). Thus, the MPs have a magnetic 
moment in this plane even in the absence of the external magnetic field. The magnetic 
moment of a platelet has the same direction as that of the total magnetic field at its site, 

m i = Ml x l y l z <k i l\k i \, (1) 

where M is the remanent magnetization of maghemite, l x , l y and l z are the dimensions of a 
platelet along the x, y and z axes respectively. With M = 50 emu/cm 3 [2j, l x = l z — 1 
and l y = 0.1 /im one obtains: rrii ~ 3.121 eV/G. 



In 



The magnetite cluster (MC) consists of nanoparticles which are 5 nm in diameter 
the case of finite temperature and finite magnetic field, the mean total moment of the MC, 
(M), is: 

(M) « = X H, (2) 

where n is the number of nanomagnets in the cluster, /x is the magnetic moment of an 
individual nanomagnet, H is the magnetic field strength at the site of the MC, T is the 
temperature and k is the Bolzmann constant. With R = 0.5 /im and r = 2.5 nm, being 
the radii of the MC and of the nanoparticle respectively, one obtains n w 8 • 10 6 . 

The proportionality constant between the magnetic moment and the field strength, x, is 
the magnetic susceptibility. Dividing it by the volume of the MC one obtains the volume 
susceptibility of the MC, Xv, which at 300 K is equal to 0.12 CGS units. With H = 10 Oe, 
which is the typical value of the total field at the site of the MC, one obtains: (JVC) ~ 0.392 
eV/G. Note that this value is about an order of magnitude smaller than the magnetic 
moment of a single MP. 

The potential energy of the MC reads as: 



E{R) = -Xv~kR 3 o 



B 



N 



(3) 



where R defines the position of the MC, B is the induction vector of the external magnetic 
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field, \imed ~ 1 is the permeability of the medium, N is the number of MPs, Hi{R) is the 
magnetic field created by the i-th MP at the site of the MC, which is known to be [1] 



H S (R) 
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(4) 



Here fj describes the position of the j-th platelet and rhj is its magnetic moment defined in 
Eq. (pQ). The total magnetic field at the site of the i— th MP is: 



B 



N 



lurried 



(5) 



The first term describes the external magnetic field while the second term describes the 
magnetic field created by all MPs except the i— th one. 

It follows from Eq. ([5]) that the total magnetic field CKj is determined by the magnetic 
moments of the platelets. Thus Eqs. (CD) and ([5]) have to be treated iteratively. In the 
zeroth-order of approximation rhi are aligned along the x-axis, what is energetically the 
most favorable configuration of the system. The magnetic moment of a platelet is then 



m 



(0) 



Ml x Ll z i, where i is the unit vector along the x-axis. The total magnetic field in the 



first-order approximation at the site of the i-th MP reads as: 



%f ] =b + wijyhta, (6) 

where is the x-coordinate of the i-th platelet and & = — Xj\ 3 . Substituting Eq. flSJ) 

into Eq. ([!]) one yields the first-order approximation for m,: 



Ml x l y l z {B x + 2Ml x l y Ui 



m\ = = (7) 

<J(B X + 2MlJ y lSf + Bl + Bl 

m« = (8) 

(1) MlJyl Z B Z 

m\ ' = ; . (9) 

(B x + 2MlJ y lS) 2 + Bl + Bl 

Here B x and B z are the x- and z-components of the external magnetic induction vector 
respectively. In [l] we demonstrated that the first-order approximation can be used to 
calculate the interaction energy with accuracy higher than 1%. 
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FIG. 2: Potential energy surfaces of the MC calculated as a function of x and y coordinates, while 
z=0 /im (see Fig.[T]) at different orientations of the external magnetic field vector: plot a): magnetic 
field vector directed along the x-axis; plot b): magnetic field vector directed along the y-axis; plot 
c): magnetic field vector directed along the z-axis. The MPs are shown with black rectangles. The 
gray rectangle in the center of the potential energy surfaces shows the region, where the MC can 
not be placed, due to its finite size. 

Figure [2] shows the potential energy surfaces of the MC as a function of its coordinates x 
and y, while z=0 \xm (see Fig.[T]), calculated at different orientations of the external magnetic 
field vector. Because of the MC size and because of the MPs there exists a forbidden region 
on the potential energy surface, where the MC can not be placed. The MPs are shown 
in Fig. [2] with black rectangles. The gray rectangle in the center of the potential energy 
surfaces defines the forbidden region for the MC. 

In our calculations the external magnetic field strength is 0.5 G, being a typical value of 
the Earth magnetic field strength. The potential energy surfaces (Fig. [2]) were calculated 
using Eq. Q. The potential energy surfaces calculated for the external magnetic field 
directed along the x-, y- and z-axes are shown in plots a), b) and c) of Fig.[2]respectively. The 
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three potential energy surfaces are similar, although some differences can be observed. The 
potential energy surfaces corresponding to the x- and z- orientations of the external field have 
axial symmetry along the y=0, z=0 axis, while the potential energy surface corresponding 
to the orientation of the external field along the y-axis has point symmetry with respect 
to the point (5.45,0) /im. There are two minima with energies about -8.5 eV at the tips 
of the maghemite chain. These minima are the global energy minima, being the spots of 
energetically most favorable attachment of the MC to the chain of MPs. This fact is in 
agreement with experimental observations, where the MC was observed at the tip of the 

n 

chain [5j. 

To illustrate the effect of the external magnetic field on the magnetoreceptor system we 
have calculated the differences in forces acting on the MC due to the 90° change of the 
direction of the external magnetic field. The force differences are shown in Fig. [3j The thin 
line shows the force difference arising due to the change of external magnetic field direction 
from x to z, and the thick line shows the force difference arising due to the change of the 
external magnetic field direction from x to y. The force differences were calculated as a 
function of x-coordinate of the MC, while y=0.8 fim and z=0 fim. Fig. [3] shows that the 
force change, caused by the 90° change of the direction of the external field is 0.1-0.2 pN in 
both cases. 

It was experimentally demonstrated j^, Q that the MCs are connected to the nerve cell 
membrane. Depending on the magnetic field strength the magnetite cluster exerts forces 
on the membrane and activates mechanosensitive ion channels increasing the flux of ions 
into the cell. The ions change the membrane potential. If the potential is reduced to the 
threshold voltage §J, an action potential is generated in the cell, which opens up hundreds 
of voltage-gated ion channels in the membrane. During the millisecond that the channels 
remain open, thousands of ions rush into the cell 8j, producing a nerve signal to the brain. 
The mechanosensitive ion channels influence the time needed for the membrane potential to 
reach the threshold value, and thus influence the birds behavior. 

A typical example of a mechanosensitive ion channel is the transduction channel of a 

nnn 

hair cell (for review see Refs. [9j,ll0j,llJ). The opening/closing of the mechanosensitive ion 
channel is regulated by the so-called gate, which is a large biological complex (protein or 
complex of proteins) at the edge of the ion channel |9j, LLC! 11, 12]. The gate is connected to 

onnn 

an elastic element, the gating spring |9J, HQ|, LLU 1121], transmitting the force to the gate. 
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FIG. 3: Difference in force acting on the MC and arising due to the 90° change of the direction of 
the external magnetic field. The thin line shows the force difference arising due to the change of 
external magnetic field direction from x to z, and thick line shows the force difference arising due 
to the change of external magnetic field direction from x to y. The force differences were calculated 
as a function of the x-coordinate of the MC, while y=0.8 /jm and z=0 /mi. 

The ion channel has two conformations: closed and open. Because the gate swings 
through a distance A upon opening, an external force / changes the energy difference between 
open and closed states and can bias the channel to spend more time in its open state. The 
gating springs are connected to the magnetite cluster which produces an external pull on 
the gates. As follows from our calculations the magnitude of this pull is about 0.2 pN, when 
the direction of the external magnetic field is changed on 90°. The work done in gating the 



channel is 12J: AE = As — fX, where the first term represents the change of the intrinsic 
energy between the open and the closed states of the channel and the second term shows 
the work of external force required for opening the channel. A is the displacement of the 
gate. For the mechanosensitive ion channels in hair cells A ~ 4 nm {{J [l^] . The probability 
for the ion channel to be open in the presence of external force is: 



1 



p= i + «p(^)- (10) 

If no external force is applied then / = and the corresponding probability for the channel 
to be open is p . Thus, the change of channel opening probability due to the applied force 
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is: 



V 



p-p _ exp(ff) (exp(§) - 1) 



Po " exp({A) +eX p(^) 




(11) 



The value of Ae is not known. Usually [12], it is assumed that Ae = 0, but in general it is 
not because the gate can form hydrogen bonds with the membrane, which break when the 
gate is opened. Thus Ae > 0. 

Fig. H] shows the dependence of the change of channel opening probability, rj on Ae (thick 
line). From Fig. H] and from Eq. (llip it follows that the change of channel opening probability 
saturates at large values of Ae. The limiting value is rj max = exp ij^) — 1. For the given /, 
A and T: r] max = 0.21, being the maximal change of channel opening probability possible in 
the suggested mechanism. If Ae = then i] = 0.096. If Ae is positive then 77 is somewhere 
between 770 and r) max . 

Another possible transducer mechanism of the geomagnetic field is based on the elastic 
deformation of the membrane. The deformation mechanism might arise in addition to the 
gating mechanism or be an alternative to it. The work performed on membrane deformation 

n 

is given by [1]: 



where AS" is the change of the membrane surface area and 7 is the membrane surface 
tension coefficient. Substituting A instead of fX in Eq. (TTTi) one obtains the change in 
channel opening probability caused by the membrane deformation. 

Fig. H] shows the dependence of the change of channel opening probability caused by 



dyn/cm=10 -5 N/m, being the typical surface tension coefficient of a membrane [13, Il4 |. 
The maximal value of i] de ^ is f]f^[ x = 0.36, being 1.7 times greater than in the case of the 
gate-spring mechanism discussed above. If Ae = then 77^ = 0.15. Since Ae is expected 
to be positive then 0.15 < 77 * < 0.36. 

In the present paper a possible mechanism of avian orientation in a magnetic field is dis- 
cussed. It was shown that in the external magnetic field the MCs experience an attractive 
(repulsive) force leading to their displacement, which induces a primary receptor potential 




(12) 



the membrane deformation, on Ae (thin line) obtained for / = 0.2 pN and 7 = 0.01 
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FIG. 4: Change of the mechanosensitive ion channel opening probability calculated as the function 
of change of the intrinsic energy between the open and the closed states of the channel. The 
thick line corresponds to the gating-spring transducer mechanism and thin line corresponds to the 
mechanism based on the elastic deformation of the membrane. 

via mechanosensitive membrane channels leading to a certain orientation effect of a bird. 
We believe that the suggested magnetoreception mechanism is a realistic candidate for the 
magnetoreception mechanism in birds. It might also be responsible for magnetosensation 
in other animals like fishes, salamanders, bees (for review see Ref. [lj]). Unfortunately, 
lack of sufficient information about magnetic particles in these species hinders us to draw 
conclusions about their precise magnetoreception mechanism. However, we believe, that the 
magnetoreception mechanism should be universal, i.e. the same for all kinds of animals with, 
probably, minor alternations. Therefore, after more experimental data regarding the mag- 
netic particles in animals become available the present investigation can be easily extended 
to a more general description. 
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